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I . INTRODUCTION 


A critical  element  in  the  prediction  of  teleseismic 
ground  motion  from  a nuclear  explosion  is  the  reflection  of 
the  compressional  wave  from  the  free  surface  above  the  working 
point.  This  reflected  phase,  pP,  interferes  with  the  direct 
wave,  P,  and  significantly  affects  body  wave  cunplitudes  at 
teleseismic  distances. 

Most  seismological  analyses  of  the  pP  event  assume  that 
the  P reflection  coefficient  is  elastic  and  that  the  source 
region  is  transparent  to  pP.  Recent  finite  difference 
simulations  of  nuclear  cratering  events  by  J.  Trulio  and 
N.  Perl  at  Applied  Theory,  Incorporated,  have  produced  a large, 
low  frequency,  negative  pulse  following  the  direct  P wave  at 
emergent  angles  appropriate  for  teleseismic  P waves.  They 
attribute  this  pulse  to  the  preferential  expansion  of  the  source 
region  toward  the  free  surface  and  the  subsequent  closure  of 
cracks  above  the  source. 

Unfortunately  there  are  few,  if  any,  near  source  ground 
motion  records  which  can  be  used  to  verify  these  assumptions. 
Both  cost  and  containment  limitations  prohibit  obtaining 
adequate  near  source  ground  motion  data  from  nuclear  explosions 
relevant  to  the  pP  issue.  Moreover,  even  if  these  limitations 
could  be  removed,  the  variation  of  ^ situ  rock  material  prop- 
erties could  alter  ground  motion  as  much  as  variations  in  pP 
when  depth  of  burial  is  changed. 

In  view  of  the  above  limitations  we  have  turned  to 
laboratory  modeling  in  order  to  obtain  ground  motion  data  as 
a function  of  depth  of  burial  under  well  controlled  conditions. 
Numerical  and  analytic  models  were  then  used  to  explain  the 
data. 
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Two  types  of  experiments  were  carried  out.  First  we 
overburied  the  charge  in  order  to  measure  the  direct  (free- 
field)  P wave  radiated  by  the  explosion.  We  will  call  this 
class  of  experiments  free-field  since  the  seismic  radiation 
is  uncontaminated  by  the  pP  event.  In  the  second  type  of 
experiment  the  charge  was  detonated  near  a free  surface  and 
the  direct  P wave  is  contaminated  by  pP.  We  will  call  this 
class  of  experiments  underburied . It  includes  as  a subset 
depths  of  burial  which  produce  craters  at  the  free  surface. 

One  of  our  objectives  was  to  determine  the  effect  of 
depth  of  burial  on  both  body  wave  and  surface  wave  magnitudes. 
This  requires  the  capability  to  analytically  continue  the 
results  of  a close-in  numerical  simulation  to  far-field  dis- 
tances. 

Rigorous  analytic  continuation  procedures  for  an 
explosive  source  detonated  in  a halfspace  are  currently  being 
developed  at  Systems,  Science  and  Software  (sM  and  were  not 
available  for  this  research  project.  Quantitative  estimates 
of  the  far-field  seismic  radiation  from  our  simulations  must 
await  the  development  of  this  continuation  procedure. 

In  spite  of  this,  the  project  has  contributed  greatly 
to  our  understanding  of  the  seismic  wave  field  radiated  by  an 
explosion  in  a half space.  In  particular,  during  the  course 
of  this  work,  we  have: 

1.  developed  and  demonstrated  the  capability 
to  design  and  instriiment  small  scale  labo- 
ratory model  experiments  and  obtain  repeat- 
able  displacement  time  history  data  from 
free-field  and  underburied  explosive  sources, 

2 . shown  that  our  numerical  models  generate 
displacement  time  histories  which  agree 
with  the  data  obtained  from  the  laboratory 
model, 
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3.  experimentally  determined  the  difference 
in  displacement  time  history  between  free- 
field  and  underburied  sources  at  emergent 
amgles  from  the  source  which  are  particularly 
relevant  to  body  wave  magnitude, 

4.  used  numerical  and  analytic  models  to  deter- 
mine the  physical  mechanism  behind  the 
generation  of  the  broad  negative  phase  which 
follows  the  direct  P arrival. 

In  Section  II  we  present  a description  of  the  laboratory 
model  and  the  results  obtained.  Section  III  compares  the 
results  of  the  numerical  model  with  the  laboratory  data.  Sec- 
tion IV  provides  the  analysis  of  the  negative  pulse.  Section 
V gives  the  conclusions  and  recommendations. 


II.  THE  LABORATORY  MODEL 

2 . 1 INTRODUCTION 

Free-field  and  underburied  cratering  explosions  were 
modeled  in  the  laboratory  on  a small  scale  by  embedding  one- 
fourth  gram  spheres  of  the  high  explosive,  PETN,  in  concrete 
cylinders  of  120  cm  dicuneter.  The  bottom  end  of  the  cylinders 
represented  the  ground  surface  and  the  seismic  coupling  of  the 
charge  was  determined  by  measuring  the  displacement  of  a point 
on  the  top  surface  directly  above  the  charge.  The  thickness 
of  the  cylinders  varied  between  33  and  60  cm.  The  vertical 
distance  between  charges  and  sensors  was  constant  (30.5  cm)  for 
all  the  experiments.  The  distance  between  the  charges  and  the 
bottom  free  svirface  was  varied  to  simulate  different  depths  of 
burial.  The  dimensions  of  the  cylinder  were  selected  such  that 
all  surfaces  except  the  bottom  surface  were  well  beyond  the 
inelastic  region  siirrounding  the  charge.  The  top  surface,  where 
sensors  were  located,  was  always  outside  the  inelastic  region 
surrounding  the  explosive. 


2.2  EXPERIMENTAL  DESIGN 

The  actual  configuration  is  shown  in  Figure  2.1  and  the 
summary  of  tests  is  given  in  Table  2.1.  Each  cylinder  was  used 
for  three  shots  fired  separately  with  charges  spaced  in  a tri- 
angular array.  The  vertices  were  at  wide  enough  separation  to 
assure  that  the  inelastic  regions  for  separate  shots  did  not 
overlap.  The  gauges  were  also  located  at  the  vertices  of  a 
triangle  with  each  gauge  directly  above  a charge.  At  Vertex  1 
horizontal  and  vertical  displacements  were  measured.  Vertical 
displacement  and  acceleration  were  measured  at  Vertex  2 and 
only  vertical  acceleration  was  measiired  at  Vertex  3.  (Only  Test 
0 deviated  from  this  configuration  in  that  two  displacement 
gauges  and  an  accelerometer  were  fixed  above  one  charge.  This 
test  was  a check  of  the  instrumentation.) 


Vertex  1 

(-15.2,  8.7,  0) 

Vertex  2 

(15.2,  8.7,  0) 

Vertex  3 

(0,  -17.6,  0) 

VI , Vertical  Displacement  Gauge 

(-14.0,  8.7,  0) 

Hi,  Horizontal  Displacement  Gauge 

(-16.5,  8.7,  -0.5) 

Sensitive  axis  | 

V2,  Vertical  Displacement  Gauge 

(14.0,  8.7,  0) 

A2 , Vertical  Accelerometer 

(16.5,  8.7,  0) 

A3,  Vertical  Accelerometer 

(0,  -17.6,  0) 

Charge  1 

(-15.2,  8.7,  -31) 

Charge  2 

(15.2,  8.7,  -31) 

Charge  3 

(0,  -17.6,  -31) 

Cylinder  Bottom 

(X,  Y,  -33  or  -37 

-60  [Table  I]) 

Figure  2.1.  Schematic  of  the  laboratory  model.  The  diameter 
of  the  cylinder  was  122  cm. 
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TABLE 


Instrument  Check 


The  concrete  cylinder  was  made  of  a material  called 
"high-early"  concrete.  This  designation  refers  to  the  fact  that 
it  attains  a high  strength  soon  after  setting.  It  has  a den- 
sity  of  2.52  gm/cm  and  unconfined  compressive  strength  of 

c 2 

345  X 10  dynes/cm  (345  bars)  with  low  porosity  and  high 
water  saturation. 

The  design  of  the  explosive  charge,  and  in  fact  the 
scale  of  the  entire  experiment,  was  fixed  by  the  need  to  use 
a source  that  was  very  nearly  spherically  symmetric.  Such  a 
source  was  available  from  Stanford  Research  Institute  (SRI) 
in  the  configuration  shown  in  Figure  2.2.  Tests  at  SRI  have 
shown  that  these  charges  yield  a shock  wave  which  is  spheri- 
cally symmetric  to  within  two  percent  [Erlich,  1975] . Each 
charge  is  a one-fourth  gram  of  PETN  powder  packed  to  a density 
of  1 gm/cm^.  The  explosive  is  encapsulated  in  a lucite  shell 
and  can  be  detonated  by  an  explosive  bridge  wire  at  the  center. 

The  principal  sensors  were  eddy-current  displacement 
gauges  supplied  by  Kaman  Sciences,  Colorado  Springs,  Colorado. 
Each  unit's  output  varies  linearly  with  distance  between  the 
probe's  tip  and  a small  piece  of  aluminum  foil  fixed  with 
epoxy  to  the  concrete  surface.  These  transducers  have  rather 
low  output  levels,  2 mV/ym,  but  they  have  the  advantages  that 
they  do  not  touch  the  concrete  and  their  response  extends  to 
50  kHz  (3  db) . To  supplement  the  displacement  gauges,  two 
accelerometers  supplied  by  Endevco  (San  Juan  Capistrano, 
California)  were  used  also.  The  accelerometers  had  a range  of 

e 2 

+ 50,000  g's  (5  x 10  m/sec  ),  a sensitivity  of  10  uV/g  and 
resonant  frequencies  above  150  kHz.  For  these  experiments  the 
acceleration  risetimes  are  roughly  4 ys  which  corresponds  to 
80  kHz.  Therefore,  the  resonance  may  be  partially  excited  and 
the  acceleration  signals  should  be  treated  as  upper  limits  for 
actual  motion  of  the  surface. 
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2 . 3 FREE-FIELD  EXPERIMENTS 

In  order  to  gain  confidence  in  the  experimental  tech- 
nique, the  first  experiment  (Test  0)  was  designed  to  demonstrate 
repeated}! lity  of  the  measurements.  Two  displacement  gauges 
and  an  accelerometer  were  in  close  proximity  directly  edaove 
one  charge  in  order  to  see  if  the  three  independent  transducers 
would  give  similar  responses  to  the  same  seismic  motion.  The 
resulting  measurements  are  shown  in  Figure  2.3.  The  data  from 
the  accelerometer  is  integrated  twice  to  give  displacement. 

The  direct  displacement  responses  were  in  fact  quite  similar, 
showing  an  arrival  time  of  92  ys,  risetime  of  15  ys,  peak 
motion  of  5 ym  and  full  width  at  half  maximtun  of  about  20  ys. 

The  double  integration  of  acceleration  data  results  in  a 
slightly  higher  displacement,  which  is  expected  because  of  the 
partial  excitation  of  resonance  of  the  acceleration  gauges 
which  was  discussed  in  the  preceding  section.  Also,  the  double 
integration  of  the  acceleration  signal  is  subject  to  a gradually 
increasing  displacement  at  late  time. 

The  results  of  free-field  Tests  4 and  5 (Table  2.1)  are 
shown  in  Figures  2.4  and  2.5.  For  these  tests,  displacement 
gauges  were  located  directly  above  each  charge.  In  addition 
to  the  primary  vertical  displacement  pulse,  which  arrives  at 
about  95  ys,  the  records  also  show  the  arrival  of  a more 
slowly  traveling  wave  at  approximately  160  ys  and  reflection 
from  the  bottom  and  sides  of  the  cylinder  at  270  ys.  The  wave 
arriving  at  160  ys,  shown  on  Hi  in  Figure  2.4,  may  be  a shear 
wave  generated  by  the  explosive.  If  this  is  the  case,  then  the 
source  for  Test  4 may  not  be  as  pure  a center  of  dilatation  as 
SRI  reported.  However,  a more  likely  explanation  is  that  the 
S waves  are  generated  by  asymmetries  in  the  emplacement  cavity 
(the  charges  were  placed  in  drill  holes  which  were  then  stemmed 
with  grout) . 

In  order  to  obtain  the  best  measure  of  free  surface 
displacement  for  these  tests  and  to  determine  the  experimental 
variability,  all  direct  measures  of  free  surface  displacement 
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MICRONS 


REFLECTIONS  FROM 
THE  BOTTOM  AND 
SIDES  OF  THE  i 
CYLINDER  7 


V 


DIRECT  SHEAR  WAVE 


• ff  (Accelerometer  #2)dt 
First  Arrival  133  ys 
+ 2.6,  - 7.6  kG's  peak 


X If  (Accelerometer  #3)dt 
First  Arrival  131  ys 
5.6,  - 8.5  kG's  peak 


Figure  2.4. 


Displacement  data  for  Test  4 where  the  charge  was 
emplaced  under  Vertex  1.  The  vertical  and  hori- 
zontal displacement  records  at  Vertex  1 and  the 
vertical  displacement  record  at  Vertex  2 are  shown 
Also  plotted  are  the  doubly  integrated  accelerom- 
eter records  from  Vertices  2 and  3.  For  the  accel 
erometer  records  the  time  of  first  arrival  and 
first  peak  and  trough  accelerations  are  indicated. 


(OP* 
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Figure  2.5.  Displacement  data  for  Test  5.  Similar  to 
Figure  2.4  except  that  now  the  charge  is 
located  under  Vertex  2. 
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from  contained  shots  were  superimposed  and  are  shown  in  Figure 
2.6.  The  repeatability  in  the  shape  and  magnitude  is  very 
good  for  these  four  records,  two  of  which  were  measured  from 
the  same  test  with  different  gauges  (Test  0)  and  two  of  which 
were  measured  from  Tests  4 and  5.  Our  estimate  of  the  total 
eunplitude  error  due  to  the  data  collection  process  was  +0.3  urn. 
The  fact  that  some  amplitude  variations  in  Figure  2.6  are  larger 
than  that  is  probedsly  due  to  slight  variations  in  the  coupling 
between  the  concrete  and  charge  surfaces.  The  results  from  Test 
5 appear  anomalous  in  terms  of  both  the  maximum  and  minimum 
amplitudes  of  the  displacement  pulse.  If  we  ignore  the  results 
from  this  test,  then  we  obtain  average  estimates  for  the  peak, 
minimum  and  static  displacements  of  4.8  ym,  0.5  ym  and  1.0  ym 
for  free-field  source. 


2.4  UNDERBURIED  SOURCES 

The  results  shown  in  Figure  2.6  give  displacement  records 
for  the  free-field  shots,  i.e.,  the  displacement  reaches  a static 
value  before  the  arrival  of  a free  surface  reflection.  We  now 
present  results  for  the  underburied  sources  where  the  reflections 
from  the  bottom  free  siirface  contaminates  the  direct  wave  from 
the  source. 

Figures  2.7  and  2.8  show  displacement  records  from  Tests 
1 and  2 (Table  2.1).  In  these  tests  the  bottom  free  surface 
was  6.5  cm  below  the  charge.  The  reflection  from  this  surface 
arrives  at  135  ym  and  is  clearly  evident  on  record  VI  in  Figure 
2.7  and  record  V2  in  Figure  2.8.  Following  the  reflected  wave, 
which  appears  to  be  an  inverted  replica  of  the  direct  wave, 
there  is  a broad  negative  pulse.  This  negative  pulse  attains 
its  minimum  value  at  210  ys  for  both  tests. 

Figures  2.9  and  2.10  show  displacement  records  from 
Tests  7 and  8 (Table  2.1} . For  these  tests  the  bottom  free 
surface  was  1.75  cm  and  1.43  cm  below  the  charge.  Craters 
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Figure  2.7. 


Displacement  data  for  Test  1.  Vertex  1 is 
directly  above  the  shot.  The  vertical  and 
horizontal  displacement  records  at  Vertex  1 
and  vertical  displacement  record  at  Vertex  2 
are  shown.  Also  plotted  is  the  doubly  inte- 
grated accelerometer  record  from  Vertex  2. 
For  the  accelerometer  records  at  Vertices  2 
and  3 the  first  arrival  and  first  peak  and 
trough  accelerations  are  indicated. 
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First  arrival  0 132  ys 
Zero  Cross  0 141  ys 


First  arrival  0 116  ys 
Zero  Cross  0 133  ys 


+12  kG's  peak  0 138  ys 
-18  kG's  peak  0 149  ys 
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- 8 kG's  peak  0 138  ys 


Figure  2.9. 


Displacement  data  for  Test  7 with  charge  located 
under  Vertex  1.  The  vertical  and  horizontal  dis 
placement  records  at  Vertex  1 and  the  vertical 
displacement  record  at  Vertex  2 are  shown.  In- 
stead of  complete  accelerometer  records  from 
Vertices  2 and  3,  times  are  given  for  the  first 
arrival  and  the  zero  crossing.  Also  indicated 
are  the  magnitudes  and  times  of  the  first  peak 
and  trough. 
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Figure  2.10.  Displacement  data  for  Test  8.  Similar  to 

Figure  2.9  except  here  the  charge  was  under 
Vertex  2. 


were  produced  at  this  surface  for  both  tests.  The  crater 
dimensions  are  given  in  Table  2.1. 


At  these  shallow  depths  of  burial,  the  broad  negative 
pulse  dominates  the  late  time  portion  of  the  seismogram.  This 
is  illustrated  further  in  Figure  2.11,  which  compares  displace- 
ment records  at  stations  directly  above  the  charge  for  Tests 
4 and  8. 


2 . 5 SUMMARY 

1.  We  have  demonstrated  that  laboratory  model  experi- 
ments can  produce  extremely  useful  and  repeatable  displacement 
time  histories  that  are  relevant  to  teleseismic  coupling. 

2.  We  have  shown  that  the  broad  negative  pulse,  which 
follows  the  P arrival  for  the  underburied  shots,  is  not  peculiar 
to  the  cratering  mode.  This  event  occurred  for  all  londerburied 
shots  including  those  that  did  not  spall  at  the  free  surface 
and/or  crater. 

3.  Our  experimental  results  indicate  that  a necessary 
condition  for  the  occurrence  of  the  negative  phase  is  that  the 
shot  be  detonated  near  the  surface  in  a half space.  Nonlinear 
material  behavior  at  the  free  surface  is  not  required. 
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III.  THE  NUMERICAL  MODEL 


3 . 1 INTRODUCTION 

Numerical  simulations  of  the  free  field  emd  underburied 
shots  were  performed  with  Lagreuigian  fin*ite  difference  codes 
using  constitutive  models  for  shear  and  tension  failure  reported 
by  Cherry,  et  , (1975),  However,  in  order  to  permit  pP  to 
arrive  as  a distinct  phase  a modification  to  the  tension  failure 
portion  of  the  model  was  required.  Therefore  we  will  present 
those  features  of  the  model  required  to  understand  the  modifica- 
tion. 

A fundamental  material  property  in  the  model  is  the  frac- 
ture locus 

Y = y (F)  = y + ^ P / 2 - , 

' n* ' 

in 

where 


Y = (3  » 


P = P 


- 1 (^') 


1/3 


J^'  and  J^'  are  the  second  and  third  deviatoric  invariants  and 
P is  the  hydrodynamic  component  of  the  stress  tensor. 

Cherry,  at  , (1975)  assumed  that  at  any  stress  state, 
F,  tension  failure  would  occur  if  the  following  conditions  were 
simultaneously  satisfied. 

1.  The  fracture  locus  was  exceeded  by  the 
unadjusted  stress  deviators. 

2.  Any  principal  stress  was  tensile. 
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This  model  produced  an  extraordinary  match  to  the  free- 
field  data.  However  pP  was  unable  to  propagate  through  the 
region  around  the  cavity  containing  open  tension  fractures. 
This  was  not  consistent  with  the  experimental  data  e.g.,  the 
result  shown  in  Figure  2.7.  Therefore  we  added  a third  condi- 
tion which  required  that  P should  be  tensile,  i.e.. 


F = < 0 


where  F is  negative  in  tension. 


(3.1) 


With  this  third  condition  the  match  to  the  free-field 
data  deteriorated  somewhat.  However  the  overall  agreement  with 
both  the  free-field  and  underburied  shots  was  more  than  satis- 
factory. 


3.2  SIMULATION  OF  THE  EXPERIMENTS 

The  calculated  ground  motion  30.5  cm  from  the  free-field 
source  is  shown  in  Figure  3.1  along  with  all  the  free-field 
experimental  data.  Two  free-field  calculations  are  presented 
in  the  figure,  one  in  which  tension  failure  was  permitted  only 
if  F was  tensile  and  one  in  which  tension  failure  could  occur 
for  any  value  of  P;  i.e. , with  and  without  Equation  (3.1) 
included  in  the  tension  failure  model.  The  reduced  velocity 
potential  (RVP)  transforms  for  the  two  calculations  are  shown 
in  Figure  3.2  and  compared  to  the  RVP  from  Test  4. 

We  see  from  these  two  figures  that  agreement  between 
simulation  and  data  is  excellent  when  Equation  (3.1)  is  not 
included  in  the  tension  failure  model.  We  then  used  this 
model  to  simulate  Test  No.  1,  with  the  charge  buried  6.5  cm 
below  the  free  surface.  The  result  is  shown  in  Figure  3.3 
which  gives  the  calculated  displacement  history  30.5  cm  below 
the  charge.  Comparing  this  calculation  with  the  data  shown 
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Figure  3.1.  Comparison  of  measured  and  numerically  simulated 
displacements . 
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. Simulation  of  Test  1 where  the  displacement  is 
calculated  30.5  cm  directly  below  the  shot. 
Note  the  absence  of  pP  which  arrives  in  the 
experimental  data  (Figure  2.7)  at  150  ys.  The 
source  region  is  opaque  to  this  event  because 
of  open  tension  fractures. 


in  Figure  2.7,  we  note  the  absence  of  the  pP  arrival  at  150  ys 
in  the  calculation.  The  negative  phase  is  present  but  not  pP. 

The  absence  of  pP  in  the  calculation  is  due  to  tension 
fractures  which  surroxmd  the  cavity  and  which  remain  open  at 
late  time.  The  late  time  crack  porosity  distribution  is  shown 
in  Figure  3.4.  This  region  is  opaque  to  pP  and  will  remain 
opaque  as  long  as  the  fractures  remain  open. 

Figures  3.5  and  3.6  compare  simulations  and  data  for 
Test  1 and  Test  8 respectively  with  Equation  (3.1)  included 
in  the  tension  fracture  model.  The  pP  event  is  now  a distinct 
arrival  in  Test  1 and  the  overall  agreement  between  the 
simulation  results  and  the  data  is  quite  acceptable. 


3 . 3 MATERIAL  PROPERTIES 

The  pressure  component  of  the  stress  tensor  for  the 
concrete  was  calculated  using 

P * Ay  + By^ 

where  y * ^ - 1.  No  ambient  air-filled  porosity  was  assumed 
. . ”o  . 

to  exist  in  the  concrete. 

The  parameters  defining  the  material  properties  for  the 
concrete  are  listed  in  Table  3.1.  These  parameters  were  derived 
from  measurements  of  P-wave  velocity,  V , density,  p,  and  the 
unconfined  compressive  strength,  Y^.  The  values  for  shear  modulus 
G,  and  bulk  modulus.  A,  were  determined  from  and  p by  assximing 
a Poisson's  ratio  of  0.1.  Additional  characteristics  were 
obtained  from  data  available  for  DF-5a  grout  (Cherry  and  Petersen 
1970)  . 

The  results  of  Finger,  et  al . , [1976],  provide  the  values 
for  the  constants  in  the  JWL  equation  of  state  for  the  products 
of  detonation  of  PETN  and  for  their  isentropic  expansion.  These 
equations  are 
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Flgvire  3.5.  Comparison  of  measured  and  predicted  displacements 
for  Test  1.  The  depth  of  burial  was  6.5  cm.  The 
gauge  was  located  30.5  cm  below  the  charge. 


FREE  SURFACE  DISPLACEMENT  (ym) 


Figure  3.6. 


Comparison  of  measured  and  predicted  displacements 
from  the  cratering  shot  (Test  8) . The  depth  of 
burial  was  1.75  cm  and  the  crater  radius  was 
approximately  8.5  cm.  The  gauge  was  located 
30.5  cm  below  the  charge. 
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,v  -RiV  . -RoV  „ 

P * A"  (1  - e ^ + B'  (1  - e 2 , (3.2) 


-R,V  -R,V 

P,  - A'e  ^ H.  B'e  . , 


(3.3) 


where  P is  the  pressure  during  the  isentropic  expansion.  For 

s 

Equations  (3.1)  and  (3.2),  V is  the  specific  volume,  1/p,  and 
the  constants  K' , B",  C",  Rj^,  R2,  and  u are  given  in  Table  3.2. 

For  the  numerical  simulation,  the  source  was  modeled  by 
a pressure  loading  function  defined  by  Equation  (3.2)  and  applied 
at  the  cavity  boundary.  The  pressure  function  was  started  with 
the  PETN  detonation  products  expanded  to  the  actual  volume  of 
the  loose-packed  PETN  which  corresponds  to  an  initial  radius, 

Rq,  of  0.39  cm.  Thus,  the  initial  pressure  was  31.2  kbar,  which 
is  well  below  the  Chapman-Jouguet  pressure  of  320  kbars  for 
normal  density  PETN. 


3 . 4 SUMMARY 

1.  A complete  set  of  calculations  have  been  performed 
with  depth  of  burial  being  the  only  variable.  The  results  of 
the  calculations  are  in  good  agreement  with  corresponding 
laboratory  data. 

2.  Enough  data  has  been  saved  from  these  calculations 
to  analytically  continue  these  numerical  results  to  the  far- 
field.  A rigorous  continuation  procedure  for  a source  in  a 
half space  is  currently  under  development  at  S ^ . 

3.  Modification  of  our  tension  fracture  model  was 
required  in  order  to  permit  pP  to  travel  through  the  nonlinear 
region  surrounding  the  cavity. 
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TABLE  3.2 


JWL  CONSTANTS  FOR  PETN 


= 7960  kbar 
= 192  kbar 


^normal 


= 6.65  kbar 
= 8.50  g/cm^ 

= 2.12  g/cm^ 

= 0.25 
=320  kbar 

3 

= 2.40  g/cm 

=1.77  g/cm^ 

= 1.0  g/cm^ 

(loosely  packed) 


r 


IV.  AN  EXPLANATION  FOR  THE  LARGE  NEGATIVE  PULSE 


In  Sections  II  and  III  laboratory  data  and  numerical 
simulations  of  explosions  in  a half space  were  presented.  All 
underburied  shots  had  two  events  in  common: 

1.  the  direct  P wave, 

2.  a large,  low  frequency  arrival  opposite 
(negative)  in  phase  to  the  direct  P wave. 

The  pP  phase  was  depth  of  burial  dependent.  A summary  of  the 
experimental  data  is  shown  in  Figure  4.1.  The  purpose  of  this 
section  is  to  present  an  explanation  for  the  negative  pulse 
and  excimine  its  effect  on  teleseismic  coupling. 

The  long  period  negative  motion  can  be  explained  as  an 
elastic  effect  due  to  near-field  reflections  from  the  free 
s\u:face.  The  negative  motion  thus  does  not  propagate  to 
teleseismic  distances.  To  demonstrate  this,  we  did  the 
following  experiment.  Using  the  reduced  velocity  potential 
from  the  free-field  source  (Figure  3.2)  as  a source  time  func- 
tion, we  propagated  the  seismic  waves  as  though  the  medium 
were  linearly  elastic.  The  source-receiver  geometry  is  shown 
in  Figure  <..2. 

The  elastic  propagation  was  done  both  numerically  and 
euialytically  to  eliminate  the  possibility  of  program  errors. 

The  numerical  calculations  were  done  with  a finite  element 
code.  The  analytic  program  is  a generalized  ray  program  which 
is  capedsle  of  propagating  the  entire  (near-field  and  far-field) 
seismic  wave  field  (Barker,  1976) . 

Figure  4.3  shows  the  comparison  between  the  vertical 
ground  motion  from  the  two-dimensional  calculation  described 
in  Section  III  with  that  obtained  by  propagating  the  free- 
field  reduced  velocity  potential  elastically.  The  elastic 
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Figure  4.1. 


The  displacement-time  histories  are  plotted  for  the 
five  experiments  to  be  further  analyzed.  The  test 
identifier  from  Table  3.1  is  noted  on  each  plot.  At 
the  right  are  the  amplitudes  of  the  peak  displace- 
ment in  microns.  On  the  plots  for  the  underburied 
experiments  we  show  the  expected  arrival  times  (from 
elastic  theory)  of  the  free  surface  pP  phase  and  the 
unwanted  reflections  from  the  sides  of  the  cylinder. 
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Figure  4.3.  The  vertical  displacement  from  the  complete 
two-dimensional  calculation  (dashed  line)  is 
compared  with  the  displacement  obtained  by 
elastically  propagating  the  free-field 
reduced  velocity  potential  (solid  line) . The 
elastic  propagation  was  performed  numerically 
using  the  RVP  of  Figure  3.2  as  a source. 


propagation  was  done  numerically  using  a finite  element  code. 

In  this  case  9-0  and  the  receiver  is  located  on  a bottom 
free  surface  at  a depth  of  30  cm.  We  see  that  the  signals  are 
very  much  alike.  In  particular,  they  both  exhibit  the  long 
period  negative  pulse  mentioned  above.  We  may  therefore  con- 
clude that  this  pulse  is  due  to  elastic  propagation  effects. 

In  order  to  understand  the  elastic  propagation  effects, 
we  turn  to  the  generalized  ray  calculations.  These  calculations 
allow  one  to  examine  separately  the  contributions  of  the  direct 
P amd  reflected  pP  and  pS  phases.  Figure  4.4  shows  sketches  of 
the  Green's  functions  (6-function  responses)  for  vertical  dis- 
placement for  a typical  source-receiver  geometry.  The  figure 
shows  the  signals  separately  and  combined.  The  direct  P wave 
begins  as  a pulse  followed  by  a plateau,  which  is  the  near- 
field P contribution.  The  pP  phase  begins  as  a pulse,  but  we 
note  that  the  near-field  part  of  this  signal  is  a negative 
ramp,  which  is  typical  of  near-field  reflections.  This  is  an 
effect  due  to  a spherical  wave  reflecting  from  a flat  impedance 
contrast,  in  this  case  the  free  surface.  The  near-field  ramp 
of  the  pS  phase  adds  to  the  pP  ramp  to  produce  a static  offset. 
For  0 » 0,  there  is  no  far-field  pS  wave,  but  there  is  a 
near-field  term. 

The  Green's  functions  were  convolved  with  the  free-field 
reduced  velocity  potential.  An  example  is  in  Figure  4.5  shown 
for  0 * 4®,  R » 30  cm.  The  shape  of  the  signal  is  very  much 
like  those  in  Figure  2.2.  Here  we  can  identify  the  long  period 
negative  pulse  with  the  negative  ramp  between  the  pP  and  sP 
arrival  times. 

To  show  how  the  waveform  develops  and  how  it  propagates 
to  teleseismic  distances , we  have  computed  the  signal  along  a 
teleseismic  ray  path,  9 = 20®,  for  increasing  values  of  R 
(Figure  4.6) . For  small  values  of  R,  the  negative  pulse  is 
not  apparent  because  the  sP  - pP  time  is  too  short.  As  this 
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Figure  4.5.  The  vertical  displacement  is  shown  for 

0 = 4®,  R * 30  cm.  The  reduced  velocity 
potential  is  propagated  using  the 
analytic  generalized  ray  program. 
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time  difference  increases,  the  long  period  negative  pulse  1 

I becomes  prominent.  At  larger  ranges,  the  pP  near-field  term 

[ becomes  small.  At  teleseismic  ranges,  the  negative  pulse  is 

I edasent  and  has  no  effect  on  measurements  of  m^^  or  M^.  These 

results  make  clear  the  importance  of  using  the  appropriate 

f < 

Green's  function  when  extrapoling  near-source  numerical 
solutions  to  the  far-field. 


V.  SUMMARY  AND  RECOMMENDATIONS 


1.  The  free-field  RVP  spectrum  for  an  explosion  in 
concrete  was  observed  to  be  mildly  peaked,  with  the  difference 
between  the  spectral  peak  and  4*^  corresponding  to  0.09  magni- 
tude units. 

2.  We  were  able  to  obtain  an  excellent  match  to  the 
free-field  data  by  using  published  data  for  the  material  prop- 
erties of  concrete,  our  own  measurements  of  P velocity  and 
unconfined  compressive  strength  and  published  data  for  PETN. 

3.  The  large,  long  period  negative  pulse  was  observed 
on  all  underburied  shots  in  both  the  laboratory  data  and  the 
numerical  simulation,  whether  or  not  spall  and/or  cratering 
occurred  at  the  free  surface.  Our  explanation  of  this  pulse 
is  that  it  is  a near-field  wave  produced  by  a spherical  P 
wave  reflecting  from  a plane  impedance  contrast.  It  should 
not  propagate  to  the  far-field. 

4 . We  were  able  to  obtain  an  acceptedsle  match  to  the 
underburied  shot  data  only  after  modifying  the  tension  failure 
model.  This  modification  somewhat  degraded  the  match  to  the 
free-field  data. 

5.  The  calculational  results  indicate  that  tension 
failure  is  required  to  produce  a peaked  spectrum.  However 
the  distinct  pP  arrival  observed  in  the  data  suggest  that 
these  fractures  either  do  not  remain  open  or  else  fill  with 
water.  It  would  be  interesting  to  section  the  concrete 
blocks  and  observe  the  degree  of  fractures  produced  in  each 
experiment. 

6.  Sufficient  data  from  each  calculation  has  been 
saved  in  order  to  analytically  continue  the  close  in  ground 
motion  to  the  far-field.  Rigorous  analytic  continuation 


algorithms  for  sources  in  a halfspace  are  being  developed  at 
S^.  They  will  be  used  to  determine  the  effect  of  depth  of 
burial  on  teleseismic  ground  motion  using  the  data  obtained 
under  this  contract. 

7.  The  laboratory  modeling  work  should  be  extended 
to  include  multiple  charge,  decoupling  and  asymmetrical 
cavity  studies. 

8.  We  should  perform  free-field  experiments  in  real 
rocks,  e.g.,  granite,  sandstone,  shale  and  salt,  both  saturated 
2und  dry  where  possible,  using  the  experimental  techniques 
developed  under  this  contract. 
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